Control of the chaotic velocity dispersion of a cold electron 
beam interacting with electrostatic waves 
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Abstract. In this article we present an application of a method of control of Hamiltonian systems to the chaotic velocity 
diffusion of a cold electron beam interacting with electrostatic waves. We numerically show the efficiency and robustness of 
the additional small control term in restoring kinetic coherence of the injected electron beam. 



INTRODUCTION 

The consequences of chaotic dynamics can be harmful in 
several contexts. During the last decade or so, much at- 
tention has been paid to the problem of chaos control. 
Controlling chaos means that one aims at reducing or 
suppressing chaos by mean of a small perturbation so 
that the original structure of the system under investiga- 
tion is kept practically unaltered, while its behavior can 
be substantially altered. In particular, in many physical 
devices there are undesirable effects due to transport phe- 
nomena that can be attributed to a chaotic dynamics. For 
example, chaos in beams of particle accelerators leads 
to a weakening of the beam luminosity [1, 2]. Similar 
problems are encountered in free electron lasers [3, 4]. 
In magnetically confined fusion plasmas, the so called 
anomalous transport, which has its microscopic origin 
in a chaotic transport of charged particles, represents a 
challenge to the attainment of high performance in fusion 
devices [5], The theoretical description of the physics of 
these devices (beam accelerators, free electron lasers, fu- 
sion devices) is based on Hamiltonian models, thus the 
above mentioned possibility of harmful consequences of 
chaos is related with Hamiltonian chaos. One way to con- 
trol transport would be that of reducing or suppressing 
chaos. There exist numerous attempts to cope with this 
problem [6]. However, these efforts have mainly been fo- 
cused on dissipative systems. For Hamiltonian systems, 
the absence of phase space attractors has been a hin- 



1 Unite Mixte de Recherche (UMR 6207) du CNRS, et des universites 
Aix-Marseille I, Aix-Marseille II et du Sud Toulon- Var. Laboratoire 
affilie a la FRUMAM (FR 2291). 



drance to the development of efficient methods of con- 
trol. Conventional approaches aim at targeting individual 
trajectories. However, when a microscopic description of 
a physical system is required, targeting of individual mi- 
croscopic trajectories is hopeless. 

The strategy we developed is based on a small but 
suitable modification of the original Hamiltonian system 
such that the controlled one has a more regular dynamics. 
In this article this strategy is illustrated through a par- 
ticular example: how to restore kinetic coherence of a 
cold electron beam interacting with electrostatic waves. 
Here chaos deteriorates the velocity "monochromaticity" 
of the injected beam and one has to find a strategy of 
control such that this effect is reduced or suppressed. 

The Hamiltonian which models the dynamics of 
charged test particles moving in electrostatic waves is 



H(p,x,t) 



Y,£ n cos(k„x-co n t + (p„), (1) 



n=\ 



where p is the test charged particle momentum and e„, 
k n , 0)„ and <p„ are respectively the amplitudes, wave num- 
bers, frequencies and phases of the electrostatic waves. 
Generically, for N > 2, the dynamics of the particles gov- 
erned by this equation is a mixture of regular and chaotic 
behaviours mainly depending on the amplitude of the 
waves. A Poincare section of the dynamics (a strobo- 
scopic plot of selected trajectories) is depicted on Fig. 1 
for two waves (N = 2). A large zone of chaotic behaviour 
occurs in between the primary resonances (nested regular 
structures). If one considers a beam of initially monoki- 
netic particles, this zone is associated with a large spread 
of the velocities after some time (left hand side panel of 
Fig. 1) since the particles are moving in the chaotic sea 




FIGURE 1. Poincare surface of section and (normalized) probability distribution functions F(p) of the velocity for the 
Hamiltonian given by Eq. (1) without control term (left panel) and with the control term given by Eq. (12) (right panel). 
The depicted numerical results are obtained in the case of two electrostatic waves using (e\,ki,0)i,tpi) = (0.045,1,0,0) and 
(32,4*2, fife, 00 = (0.045,1,1,0). Reprinted from Ref. [11]. 



created by the overlap of the two resonances [7]. Here 
the problem of control is to find a control term / given 
by 

9 

f(x,t) = £ J7„cos(fc n x-£V + *«), (2) 
«=i 

that is the right set of additional controlling waves, such 
that the controlled Hamiltonian 

H c (p,x,t) =H(j>,x,t)+f(x,t), 

given by the original one plus the control term, is more 
regular than H. In general, adding a generic perturbation 
adds more resonances and hence more chaos. Neverthe- 
less this problem has some obvious solutions with addi- 
tional waves having amplitudes of the same order as the 
initial ones. For energetic purposes, the additional waves 
must have amplitudes which are much smaller than the 
initial ones (max„ |tj„| <C max„ |e„|). This means finding 
a specifically designed small control term that is a slight 
modification of the system but which drastically changes 
the dynamical behaviour from chaotic to regular. In par- 
ticular one aims at building barriers to transport that pre- 
vent large scale chaos to occur in the system. 



LOCAL CONTROL METHOD AND 
CONTROL TERM 

In this section we briefly sketch the local control method 
that was extensively discussed in Refs. [8, 9]. Let us 



consider Hamiltonian systems of the form 

H(A,d)=H (A) + eV(A,d), 

that after a suitable expansion can be rewritten as 

ff(A,8) = fi)-A + ev(0)+w(A,8), (3) 

where (A, 8) 6 M 1 x T L are action-angle variables for 
Ho in a phase space of dimension 2L. The last term of 
the expansion can be written as 

w(A,8) = ewi(fl) -A + vv 2 (A,e), 

where W2 is quadratic in the actions, i.e. W2(0, 6) = 
and o|aW2(0) 0) —0. We assume that 0) is a non-resonant 
vector of M. L , i.e. there is no non-zero k G Z such that 
CO ■ k = 0. Moreover we assume without restrictions that 

j jL v(e)d L e = o. 

We consider a region near A = 0. The controlled 
Hamiltonian we construct is given by 

H c (A,0) = (Q-A + ev(e) + w{A,e) + e 2 f(e). (4) 

The control term depends only on angle variables and its 
expression is given by 

/(e) = -£- 2 w(-er<3 e v,e), (5) 

where dgv denotes the first derivatives of v with respect 
to 6 : dgv = Lkez L z"kvke*' 8 and where the linear oper- 
ator r is a pseudo-inverse of CO ■ dg, i.e. its action on a 



function v(6) = Lkez L is : 



rv(fl) 



Vk -e !k e . 



keZ L \{0} 



id k 



(6) 



We prove (see Ref. [9]) that H c has an invariant torus 
located at A = — eFdgv. For Hamiltonian systems with 
two degrees of freedom, such an invariant torus acts as 
a barrier to diffusion. For the construction of the control 
term, we notice that we do not require that the quadratic 
part of w is small in order to have a control term of order 



Control term for a two wave model 

We consider the following Hamiltonian with two trav- 
eling waves: 



H(p,x,t) = — +£icos(kix-(0it + (pi] 
+e 2 cos(k 2 x-co 2 t + <Pi), 



(7) 



where the wavenumbers are chosen according to a given 
dispersion relation k\ = K(cO\) and k 2 = K(co 2 ). The 
numerical results we will show in the following have 
been obtained using (ei,k\,COi,cpi) = (e, 1,0,0) and 
(£ 2 ,k 2 , tBJi tyi) = ( e , 1,1,0). Similar results are expected 
for other values of (e n ,k n ,(O n ,(p n ). 

Figure 2 (left panel) depicts the probability distribu- 
tion function of the momenta of a trajectory of Hamil- 
tonian (7) for t G [0,80071:] (with initial condition cho- 
sen in the chaotic sea, e.g. for (p w 0.35,x = 0)) as a 
function of the amplitude e of the perturbation. It shows 
that after £ c w 0.028 there is no longer any barrier in 
phase space. We notice that the value of £ c for which 
the last invariant torus is broken is approximately equal 
to e c w 0.02759 [10]. 

The first step in order to compute the local control 
term which recreates a barrier in phase space is to map 
this Hamiltonian with one and a half degrees of freedom 
into an autonomous Hamiltonian with two degrees of 
freedom by considering that t mod 27T is an additional 
angle variable. We denote E its conjugate action. The 
autonomous Hamiltonian is 



H = E + — + eicos(k l x- a>it + (jh) 
+E2Cos(k 2 x - cotf + (fa). 



(8) 



Then, the momentum p is shifted by CO in order to define 
a local control in the region p«0. The Hamiltonian is 
rewritten as 



H = E + cop + £iCOs(fciJC - (d\t + (pi] 

P 2 

+e 2 cos (k 2 x - cotf + <P2 ) + ■ 



We define 

H =E + COp, 
v(x,t) = E\ cos(fcix — co\t + q>i) + s 2 cos(k 2 x— a^t + <p 2 ), 

w{p,x,t) =p 2 /2, 

provided (0 ^ (to\jk\ and CO ^ 02/^2- in this cas e the 
control term is given by 

f{x,t) = -\{Td x v)\ 

where d x is the partial derivative with respect to x and the 
action of T follows from Eq. (6). Therefore the control 
term is equal to 



f(x,t) 



£\k\ 



cos{k\x— co\t + cpi) 



+ e 2 k 2 



coki - coi 
cos(k 2 x— cxyit + cp 2 ) 



1 2 



ft)fc2-C02 



(10) 



Adding the exact control term (10) to the Hamilto- 
nian (9), an invariant KAM torus with frequency CO is 
recreated. This barrier prevents the electron beam to dif- 
fuse in phase space and the electron kinetic coherence is 
restored. In Fig. 2 (right panel) the recreation of barrier in 
phase space can be observed on the probability distribu- 
tion function of the momenta of a trajectory of Hamilto- 
nian (7) plus the control term (10) for t G [0, 8007T] (with 
initial condition chosen as in the uncontrolled case, e.g. 
for (p :=y 0.35,x = 0)). In fact for all the values of the am- 
plitude e of the perturbation, the kinetic coherence of the 
initial electron beam is preserved. 

The local control ensures that the barrier persists for 
all the magnitudes of the perturbation and also gives us 
explicitly the equation of the recreated KAM torus, that 
is 



Po(x,t) = CO-Eih 



COS^X — C0\t + cpi) 

cok\ - co\ 



, cos{k 2 x - oty + cp 2 ) 
-e 2 k 2 — i — 212. (11) 

ft>£2-®2 

The control term (10) has four Fourier modes, 

(2fci,-2©i), (2k 2 ,-2an), ((*i + k 2 ),-(an + an)) 

and ((k\ — k 2 ),—{C0\ —02))- If we consider CO close 
to {co\/k\ + O2/&2 )/2. the main Fourier mode of the 
control term is 



h - 



£162^1^2 



2(0)ii-0)i)(fl)i2-fi^) 
COs[(fei +k 2 )x- (C0i + C02)f + <Pi + P2] -(12) 



Assuming co = (C0\/k\ + Ciy 2 /k 2 )/2, the control term 
given by Eq. (12) is 



(9) h - 



2£i£2 
(vi - V 2 ) 2 



cos[(fei + k 2 )x - (ffli + G>i)t + <pi + cp 2 }, 
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FIGURE 2. Distribution of momenta p as a function of e for Hamiltonian (7) (left panel) and for Hamiltonian (7) plus the control 
term (10) (right panel). The numerical results are obtained for the following values of the parameters: (e\,k\,a)\ , (p\) = (e, 1,0,0) 
and (£ 2 ,k 2 ,a>2,<P2) = (e, 1,1,0). 
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FIGURE 3. Distribution of momentum p as a function of £ 
for Hamiltonian (7) with the approximate control term (12). 
The numerical results are obtained for the same values of the 
parameters as in Fig. 2. 

where v„ is the group velocity of the wave n defined as 
v„ = (O n /k n . 

Figure 3 depicts the probability distribution function 
of the momenta of a trajectory of Hamiltonian (7) with 
the approximate control term (12), for t G [0,8007c] and 
with initial condition (p « 0.35,x = 0), as a function 
of the amplitude e of the perturbation. It shows that 
even after e c « 0.028 barriers in phase space have been 
created. The approximate control is efficient till £ c w 
0.07. 

Robustness of the control: Influence of a 
change of frequency and phase in the 
control term 

The main problem is that the wavenumber k\ + k 2 of 
the control term does not satisfy in general the dispersion 



relation k = K(a>), i.e., we do not have in general k\ + 
k 2 = K(coi + ©2) since the dispersion relation is not 
linear. Therefore the determination of the frequency and 
wavenumber of the control term contains errors. The 
approximate control term has the form 

Am) = SiSihk 2 x 

2 2((ok\ - (Oi)(cok 2 - ©2) 

cos[(fci + k 2 )x - (©1 + ©2 + 8co)t + <pi + (p 2 \. 

We also test the robustness of the local control when 
there is an error on the phase, i.e. we consider an ap- 
proximate control term of the form 

Jq>) £i£ 2 ^1^2 

2 2(<Qki — <Qi)(<Dk2 — Oh) 

cos[(fci + k 2 )x - (©! + ©2)f + (pi + (p 2 + 8(pj. 

Figure 4 depicts the probability distribution function of 
the momenta of a trajectory of Hamiltonian (7) with the 

approximate control term (left panel), or with 
(right panel), for t G [0, 8007c] and with initial condition 
(p s=s 0.35,x = 0), as a function of the amplitude e of 
the perturbation. From these figures we notice that it is 
important to adjust precisely the frequency of the control 
term, and also choose a set of modes such that k\ + k 2 is 
close toK(a>i + ©2). 

In summary, we have shown here how our method of 
control works in suppressing chaotic velocity diffusion 
induced in a cold electron beam interacting with electro- 
static waves. These numerical results have been success- 
fully confirmed by an experimental check performed on 
a Traveling Wave Tube [11]. According to the fact that if 
the amplitude of the potential is of order e then the am- 
plitude of the control term is of order e 2 , the control term 
is realized with an additional cost of energy that is less 
than 0.1% of the initial electrostatic energy of the two 
wave system. 
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FIGURE 4. Distribution of momentum p as a function of £ for Hamiltonian (7) with the approximate control term (12) containing 
an error of Sco = 0.001 in the frequency (left panel) and an error of Sep = 0.2 in the phase (right panel). The numerical results are 
obtained for the same values of the parameters as in Fig. 2. 
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